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1 INTRODUCTION 



ABSTRACT 

In this paper, w e hav e extended the entropy-driven model of cluster evolution devel- 
oped by Bower (1997) in order to be able to predict the evolution of galaxy clusters for 
a range of cosmological scenarios. We have applied this model to recent measurements 
of the evolution of the L x — T normalisation and X-ray luminosity function in order 
to place constraints on cosmological parameters. We find that these measurements 
alone do not select a particular cosmological frame-work. An additional constraint is 
required on the effective slope of the power spectrum to break the degeneracy that 
exists between this and the background cosmology. We have therefore included a the- 
oretical calculation of the dependence on the power spectrum, based on the cold 
dark matter paradigm, which infers £!o < 0.55 (0.1 < fio < 0.7 for f2o + Ao = 1), at the 
95% confidence level. Alternatively, an independent measurement of the slope of the 
power spectrum from galaxy clustering requires f2 < 0.6 (fio < 0.65 for fl + A = 1), 
again to 95% confidence. The rate of entropy evolution is insensitive to the values of 
Qq considered, although is sensitive to changes in the distribution of the intracluster 
medium. 

Key words: galaxies: clusters - cosmology: theory 



ray emission by some other means of mass estimation, such 



Clustery uf gaLuues die the LngesL vhialised 



tions inl Uil piLJLiiL duy UuiuLUL. TIiuj, nululiuiiui^ jliid 

■ ies offe r a unique method for directly determining the rate 
!_h ' at which such structures grow in mass. This is influenced by 
the competing effects of the steepness of the density fluctu- 
ation power spectrum, characterised by the effective slope, 
n, and the current values of the cosmological parameters, 
fio and Ao In this paper, we explore the possibility that 
the value of fio, or the combination, fio + Ao = 1 (which 
produces a flat geometry), can be robustly determined from 
an approach based solely on the X-ray evolution of galaxy 
clusters. In order to be successful however, such an approach 
needs to overcome two substantial obstacles. Firstly, the 
correspondence between cluster mass and X-ray luminos- 
ity is not direct, being sensitive to the way in which the 
gas fills the cluster's gravitational potential well. Secondly, 
the degeneracy between the cosmological parameters and 
the effective slope of the power spectrum on cluster scales, 



as the lumin osity-temperature relation (e.g. Mushotzky & 
Scharf 1997) or gravitational lensing effects ( |S mail ct al 



1997, Bower & Smail 1997). The second problem can be ap- 
proached in several ways. The most straightforward is to 
adopt a fluctuation spectrum on the grounds of a physical 
hypo thesis, for example t he cold dark matter (CDM) model 



(e.g. |Bardeen et al. 1986|) . Alternatively, an empirical mea- 



surement of the power spectrum could be used, for example 
using the large-sc ale distribution of rich clusters from red- 
shift surveys (e.g. Tadros, Efstathiou fc Dalton 1997 ) or by 
derivation from the sh ape of the cluster temperat ure distri- 
bution function (e.g. Oukbir & Blanchard 1997). Perhaps 



be broken before a unique value can be selected. 



the most appealing approach would be to use X-ray obser- 
vations spanning a wide range in redshift and luminosity to 
separate out the models purely on the basis of their observed 
evolution. 

Numerous approaches to the determination of fi p have 
already been presented in the literature (e.g. see Dekel, 
Burstein & White 1996, for a recent review), which range 



The first problem can be tackled by measuring the rate 
of luminosity evolution and calibrating the efficiency of X- 



* Ao = A/ 3Hq , with the Hubble constant taking the value Ho 
lOOhkms -1 Mpc" 1 



from direct arguments based on the peculiar motions of 
galaxies in the local Universe, through indirect methods us- 
ing the baryon fraction in clusters of galaxies, to measure- 
ments of the acoustic peak in the microwave background 
spectrum. Each of these approaches has their own strengths 
and measures a different aspect of the overall cosmological 
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model. A value based on the rate of cluster mass growth is 
appealing since it measures Q.q on the basis of its large-scale 
effect over a modest factor of the Universe's expansion. Con- 
vergence of all these methods will act as confirmation that 
our global cosmological picture is valid and that no crucial 
additional physical processes have been omitted. 

Discussion of the implications of cluster evolution is also 
timely given the growing area of sky that has now been ex- 
ploited in X-ray surveys. These range in strategy fro m wide- 
area pr ojects based on the ROSA T all-sk y survey (Ebeling 
et al. 1 ^971, |Po Grandi et al. 19971 ), EMSS ( |Henry et aT. 1992| 



Gioia fc Luppino 1994), and the ROSAT North Ecliptic Pole 
survey ( |Hcnry ct al. 1999 ), to smaller solid-angle surveys 



based on s erendipitous source s identified in deep ROSAT 
fields ( e.g 



Rosati et al. 1998 



Scharf et al. 1997, Collins 



et al. 1397). Some constraints are also available from very 
deep survey fields ( Hasinger et al. 1998, McHardy et al. 1998 
Bower et al. 199€| ) although the area covered by these fields 



is currently very small. In addition, reliable temperature 
data is also becoming available for clusters spanning a range 
of redshifts and luminosities, via the ASCA satellite (e.g. 



Tsuru et al. 199rj , |Mushotzky fc Scharf 1997| , |Markevitch 
1998]). these" recent advances have motivated our study, but 
there is also a clear need to identify the places where future 
observations, based for example on the forthcoming AXAF 
and XMM missions are best targeted. One crucial question 
is whether more progress is to be made by going to lower 
flux levels, or more uniform surveys covering a wider area of 
sky. 

Our work is related to that of several others in recent 
literature, representing a range of possible approaches to the 
problem of cluster evolution. This paper uses a phenomeno- 
logical model that separates out the evolution of clusters 
into factors depending on both the evolution of the mass 
spectrum and processes resulting in heating and cooling of 
the intracluster gas. The foundatio ns of the model were dis- 
cussed extensively in Bower 1997] (hereafter Paper I), pro- 
viding both a physical basis and interpretation. It allows us 
to minimise additional theoretical input into the calculation, 
by using simple scaling relations to translate the properties 
of the cluster sample at low redshifts into their equivalents at 
earlier epochs. The appro aches used by Mathiesen & Evrard 
( h998|) , Reichart et al. ( |l998| ) and Blanchard & Bartlett 
(1998) are related, but use an empirical model for the X- 
ray luminosity calibration and evolution, combined with the 
Press-Schechter method (Press & Schechter 1974) for the 
distri bution of cluster masses. The approach of Kitayama & 
Suto (1997) is more different in the sense that they achieve 



a match to the luminosity function data through varying 
the epoch at which the clusters must form. In the case of 
very low values of Qo, the distinction between the epoch at 
which a cluster is observed and that at which it is formed 
becomes important. We have therefore developed the model 
from Paper I to incorporate this effect. Finally, a method 



of direct deconvolution has been outlined by Henry (1997 



see also Eke et al. 1998) using recent ASCA temperature 



data. A wide variety of approaches to this topic is clearly 
desirable in order to indicate the robustness of the underly- 
ing principles. We have therefore included in our discussion, 
a comparison between our work and others and we outline 
the areas of uncertainty that can be considerably improved 
from further observations. 



The layout of this paper is as follows: in §2, we out- 
line the model on which this paper is based and show how 
it can be readily extended to incorporate evolution in both 
open (tt < 1, A = 0) and flat (Q. a < 1, tto ± A = 1), 
sub-critical density universes. We have also included in this 
section, our approach to incorporate the effects of the clus- 
ter formation epoch. §3 summarises the constraints on the 
model parameters using currently available X-ray cluster 
data for a range of cosmologies and we investigate the limits 
that can, at present, be set by this approach. We also fo- 
cus on the inherent degeneracy that exists between S7o and 
n, and investigate whether we can empirically distinguish 
between different cosmological models, using the methods 
mentioned earlier. The possibility of measuring evolution 
of the cluster core radius to place further constraints on 
the model is also discussed. In §4 we summarise our results 
and investigate the robustness of assumptions made in this 
model. We also compare our results with the values that 
have been obtained by other authors, and explore the differ- 
ences between the proposed models leading us to consider 
the overall accuracy of the method. With this in mind, we 
identify key strategies for future X-ray surveys. Finally, in 
section 85 we reiterate our conclusions. 



2 MODELLING THE X-RAY EVOLUTION IN 
DIFFERENT COSMOLOGICAL SCENARIOS 

2.1 X— ray Emission and the Cluster Core 

It is now well established that X-ray emission from the 
intracluster medium is dominated by collisional processes. 
For hot clusters, the most important of these is thermal 
bremsstrahlung. The emissivity scales as p 2 T a , where we 
use a — 0.4 [±0.1] for bolometric/ wide-band detectors and 
q = [±0.1] for low energy band-passes, although we find 
that the results are not affected when making changes to 
these values within the limits quoted in square brackets. 
Detections of galaxy clusters in this region of the electro- 
magnetic spectrum are sensitive to the way in which the in- 
tracluster gas is distributed. Surface brightness distributions 
are generally fitted using the following (so-called /3 — model) 
density profile 



p(r) 



Pc 



(3 



(2.1) 



with r c defining the effective core-size and f3 the rate at 
which the density falls off with radius. We adopt a constant 
value of j3 = 2/3 for the main results in this paper. This is 
appropriate for a non-singular isothermal distribu tion and 
is in agreemen t with the observational average (e.g. Jones & 
Forman 1984). The effect of departures from this asymptotic 



slope is discussed in §4.1. The emission is thus characterised 
by a flattening at small scales (typically r c ~ 100 kpc), de- 
parting from the distribution expected if the gas followed 
the underlying dark matter. A plausib le ph ysical interpre- 
tation was given by Evrard & Henry (1991, hereafter EH, 



see also Kaiser 1991 ) , hypothesising that the intracluster gas 
was pre-heated prior to the cluster's formation and has re- 
tained the entropy acquired from this pre-collapse phase. 
Assuming that any dissipation was negligible, this provides 



© 0000 RAS, MNRAS 000, 000-000 



Measuring Q from the Entropy Evolution of Clusters 3 



an entropy floor, s m in [I], forcing the gas to build up a mass 
distribution that reflects this constraint. For an isothermal 
distribution, s m i n directly corresponds to the core density, 
p c . In paper I, this idea was developed further by allowing 
Smin to evolve, using the parameterization 

Smin = s min (z = 0) + Cv eln(l + z), (2.2) 

where e determines the rate of core entropy evolution, which 
could be dominated by, for example, the gas being shock- 
heated as a result of merging (implying negative values) or 
conversely, radiative cooling (positive values) . The value e = 
corresponds to no net evolution of the core entropy. 



and number density, n oc (1 + z) ™+ a . 

The above relations are not appropriate for the scaling 
properties of individual clusters. The histories of particu- 
lar objects will differ wildly. Rather, we use only the Weak 
Self-Similarity Principle to apply these scalings to the mean 
evolution of the population as a whole (c.f. Paper I). This is 
strictly correct if the universe has critical density and a scale 
free spectrum of density perturbations. Note, however, that 
measurements based on a single cluster population (e.g. the 
slope of the L x — T relation) cannot be constrained using 
this approach since the formation histories of clusters with 
greater or smaller masses than the characteristic value may 
be different. 



2.2 Evolution in an Q — 1 Universe 



The case of evolution in a critical density universe was dis- 
cussed extensively in paper I. We outline the details again 
here before going on to discuss evolution in the general cos- 
mological context. The X-ray luminosity of an individual 
cluster is calculated by integrating over the virialised region 



£x oc 



1 p 2 T a dr. 



(2.3) 



We use Eq. 2.1 to describe the density profile, and assume 



that the temperature profile of the gas distribution can be 
described in a similar way, with the characteristic tempera- 
ture in proportion to the virial temperature of the system. 
In practice, the integral is dominated by the contribution 
from within a few core radii, and thus the scaling proper- 
ties of this integral depend weakly on the assumed density 
profile. Furthermore, departures from the standard profile 
can be accommodated by redefining the core radius of the 
system. As in EH, extracting the scaling properties from the 
integral produces the following result: 



L2 3 t-io 
x oc p c r c 1 



(2.4) 



Using Eq. 2.2 to parametrise the evolution of the core 



entropy, the core density evolves as 

p c OcT^ (l + s)-TT. 



(2.5) 



The core density can also be related to the virial radius, 
virial density and the core radius, using the asymptotic pro- 
file of the gas distribution, p c ~ p v i T (■^ i£ ) 3 ' 3 - Note that 
the effectiveness of a certain value of e is dependent on the 
/3-profile adopted. Coupled with the definition of the virial 
radius, rvir oc (Mp~ ir )3 (where M is the virial mass of the 
cluster) this leads to a relation for the characteristic cluster 
core radius 



r c oc Ms T 2 f* (1 + z) 



-l+7f(l+f) 



(2.6) 



We combine these results with the scaling relations for mass 

-6 2 

(M oc (1 + z) ™+ 3 ) and virial temperature (TocAfa (l + z)) 
to obtain the redshift dependence of the cluster luminosity, 



L B( x(l+*) !K *< 1+ * ) - 



(2.7) 



' we use the definition of specific entropy, s = c u ln(Tp 1 7 ) ; 7 is 
assumed to be 5/3, appropriate for a non— relativistic ideal gas 



2.3 Extension to flo < 1 Cosmologies 

For a sub-critical density Universe, several modifications 
are required. Firstly, the Weak Self-similarity Principle can- 
not be rigorously justified in a low density universe since 
the internal mass structure of the clusters need not be ho- 
mologous and will depend on the epoch at which they col- 
lapse. In iV-body simul ations, however, deviations f rom ho- 
mology are small (e.g. Eke, Navarro & Frenk 1998, for the 



fio = 0.3, Ao = 0.7 scenario); we will assume that any depar- 
ture from homology can be incorporated into the definition 
of the e parameter. For the aim of this paper, ie. compar- 
ing the evolutionary properties of clusters in critical and low 
density universes, this approach is adequate. 

The characteristic density of the cluster system will 
however depend on the epoch at which the system col- 
lapses (zf). For a very low density Universe, this may be 
significantly different from the epoch at which the cluster 
is observed. To relate the dynamical properties of cluster 
populations to the background properties set by the cos- 
mology, we use the results of the spherical top-hat collapse 
model (e.g. Lahav et al. 1991). This predicts that that over- 



dense spherically symmetric perturbations depart from the 
linear regime, turn around and collapse, forming virialised 
structures with mean internal densities given by the formula 
pvir = Avir Pcrit , where p cr it is the density required to close 
the Universe. The quantity A v i r is a function of the back- 
ground cosmology and ca n be approximated by using the 
following fitting functions (Eke, Navarro fc Frenk 199q): 



A vir = 178 n u 
= 178 O 



(Ao = 0) 

(fio + Ao = 1). 



(2.8) 



Finally, the linear growth factor, 8+, evolves more 
slowly than (1 + z) _1 f or Op < 1 mo dels. The appropriate 
relations are as follows (Peebles 1980) 



fio < 1 : - 
8+(x) 

x 

Qo + Ao = 1 : - 
5+ (y) 



, 3 3(l + :r)2 1 
1 + - + — 3-^-hi (1 + x) 2 

x x? 

u 



x? 



-1 




l + z 



(y 3 + 2)i 



u 3 + 2 



du. 
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y = 



2(O " 1 -i)i 



1 + 2 



(2.9) 



Accounting for these differences produces the following scal- 
ing relation for the X-ray luminosity: 



L x oc T a (l + z)"{l + z f ) 



-1 H 



A vir (z/) 



Q(z f ) 



(2.10) 



1 3 



d 



i _ ^ 
/3 2' 

The temperature scaling relation is thus 

i 

Toe 5+( 2 )^T3 (1 + Z/ ) ' V 



A vir (z/) 



0(2/) 

and the number density scales as S 



(2.11) 



6 
,-3 



2.4 The Epoch of Cluster Formation 

Since the linear growth of fl uctuations 'freezes out' when 



o- 



1 (Peebles 1980), producing the present-day 



abundance of clusters requires structure to have formed at 
progressively earlier epochs for lower density Universes. Con- 
sequently, the epoch of cluster formation, z/ can be signifi- 
cantly different from the redshift range of an observed sam 
pie, leading to inaccu racies in the scaling properties (Ki- 
tayama & Suto 1996). We tackle this problem by defining 
the epoch of formation to be such that the cluster has ac- 
quired a fraction, /, of it's total mass when observed 



M(z f (z)) = fM(z). 



(2.12) 



Material that has been accumulated between Zf and z is 
assumed to have negligible effects on the normalisation of the 
mass distribution, even though the total mass and therefore 
temperature of clusters can change. 

The effect that varying / has on the X-ray luminosity 
is illustrated in Figure [l] for a typical median redshift of 
current distant cluster samples (z = 0.4). This is done by 
plotting the following function 



F(O ,A ,2) = F (l + z f )- 



A vir ( 2/ ) 



Q(z f ) 



(2.13) 



varying with Oo, i.e. the part of equation 2.10 that depends 
on Zf. The constant, Fq, normalises the relation to the cor- 
responding value at the present day (zf(z — 0)). Clearly 
for the 0=1 case, the formation epoch scales in direct pro- 
portion to the observed redshift, with the normalisation only 
depending on the slope of the power spectrum, n (which also 
controls the rate at which structure grows). As the matter 
density drops below the critical value, we see a decrease in 
about a factor of 2 for Oo = 0.1, although the dependence 
on / is relatively weak (even for O = 0.1) : reducing / from 
1 — ► 0.5 produces only a 40% change in F, even in this ex- 



1 


_«L i ^ s - J — — 










Y 








1 


(a) f) <1.0 ; A Q =0 

1,1,1 



0.2 0.4 0.6 0.8 1 

0„ 




a, 



Figure 1. The X-ray luminosity factor that depends on the clus- 
ter formation redshift, Zf, as a function of the background cos- 
mology, for both open (f2o < 1) and flat (Oo + Ao = 1) models, 
for the case 2 = 0.4. The solid line illustrates the variation of F 
for a mass fraction of 1 and the dashed lines show progressively 
lower values, in steps of 0.1, to / = 0.5 . 



treme case. Therefore, we have selected a fiducial value of 
/ = 0.5 for the results that follow. 



3 USING THE EVOLUTION OF CLUSTERS 
TO CONSTRAIN O 

The essence of our approach is to scale characteristic prop- 
erties of a cluster population from one epoch to another. 
We therefore need to use data based on local samples as the 
basis for the scaling transformations described in the last 
section. We can then determine the likelihood distribution 
of parameters by fitting the scaled relations to the data avail- 
able at higher redshift. At this stage, we have assumed that 
our free parameters are the slope of the power spectrum, 
n, and the entropy evolution parameter, e. We have cho- 
sen to fix the cosmological density at four values (with and 
without the cosmological constant) : Oo = 0.1,0.3,0.5,1.0, 
in order to clearly illustrate the effects of a varying cos- 
mological background on the physical evolution of clus- 
ters. Below, we discuss the data-sets used to constrain the 
set (n, e) and the corresponding results. All data presented 
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Parameter 


Best Fit 


95% (min,max) 


n 


-0.10 


(-0.63,-1-0.43) 


A 


0.292 


(0.325,0.260) 



Table 1. Values of the parameters used to fit the (temperature) 
evolution of the L x — T relation, assuming a power law dependence 
on L x with slope, A and an evolution term between low and high 
redshift parameterised as (1 + z) v . Shown are the best— fit values 
of the slope and evolution parameters and 95% confidence levels 
on their dispersion, for the median redshift of the high-z sample, 
(z) = 0.3, and a value of qo = 0. 



Survey 

SHARC feurkc ct al~997|) 
WARPS fjoncs ct al. 199 



RDCS (Rosati ct al. 1998^ 



EMSS If Henry ct al. 1992 ) 



EMSS 2(Luppino & Gioia 199E) 



0.44 
0.47 
0.60 
0.33 
0.66 



N c i 

16 
11 
11 

23 
6 



Table 2. Details of the distant cluster survey XLF's used for con- 
straining evolution of the ROSAT BCS XLF. The EMSS 2 data 
is not binned (*) but rather taken as a point on the cumulative 
luminosity function. 



have been compiled with an assumed Hubble constant of 
H = SOkms^Mpc" 1 . 



3.1 The Luminosity— Temperature Relation 

We place a constraint on the temperature evolution of clus- 
ters by making use of the L x — T relation, which is at present, 
adequately described by a fixed power law, L x oc T A , and 
fixed intrinsic scatter. We determine the evolution of the 
L x — T relation by fitting a maximum likelihood model to 
the comb ined low and high redshift data given by David et 
al. ( 1995 ) and Mushotzky & Scharf ( 1997 ), compiled mainly 
from the ASCA and Einstein satellites. We assume the dis- 
tribution of temperatures (that lead to considerable scat- 
ter in the relation) are Gaussian distributed, and convert 
the 90% systematic errors in the temperature measurements 
(quoted by the authors) to their equivalent la values. In 
order to ensure that there was good overlap in luminosity 
between the high and low redshift data-sets, we limited the 
comparison to clusters with bolometric luminosities greater 
than 10 44 ' 5 ergs -1 , although the results are not particularly 
sensitive to this choice. The model for the L x — T correlation 
that we fit includes an adjustable zero-point, slope and in- 
trinsic scatter, as well as a redshift -dependent normalisation 
term that is parameterised as A log To = f] log(l + z) (where 
To is a reference temperature) . We determine confidence lim- 
its on the evolution of the normalisation by minimising over 
the other parameters and calculating the distribution of C 
values, where 



N 



(3.1) 



with the subscript i running over each cluster, Pi is the prob- 
ability of measuring the cluster with a given temperature, 
luminosity and normalisation, controlled by 77. We have as- 
s umed that AC is distributed as \ 2 with one free parameter 
( |Cash 1979| ). 

This treatment assumes that we are not interested in 
the slope of the relation. This is exactly true if the lumi- 
nosities of the clusters do not evolve and the evolution of 
the L x — T relation is due to the temperature evolution of 
the clusters alone. In practice, however, a particular choice 
of the parameters (n, e) implies correlated changes in lumi- 
nosity and temperature. For simplicity and consistency with 
other work, we have converted the evolution in both quan- 
tities into a change in temperature at fixed luminosity. This 
correction involves the slope of the relation, although we 



emphasise that this is an artefact of the way in which the 
normalisation is quoted. In principle, each model requires 
that we examine the likelihood as a function of a combina- 
tion of rj and the fitted slope, A (which changes with the 
choice of (n,e) parameters); in practice, however, the slope 
of the relation is sufficiently constrained that almost iden- 
tical results are obtained if we examine the likelihood as a 
function of the single parameter rj and then simply adopt 
the best fitting slope in the model calculation. 

The limiting values of the evolution rate parameter are 
given in Table [D along with the corresponding slope of the 
L x — T relation, for the median redshift of the high-redshift 
sample ((z) = 0.3) and a value of qo — (Qo = Ao = 0). 
These values are converted to the appropriate cosmologies 
when required, since the normalisation is affected by the as- 
sumed value of Qo through the distance dependence of the 
luminosity. In Qo < 1 cosmologies, the luminosities of clus- 
ters are brighter, pushing the evolution of the temperature 
normalisation in the positive direction for higher values of 
Qo. The uncertainties in the evolution rate are dominated 
by the intrinsic scatter in the relation. We therefore note 
that while the evolution of the relation is statistically well 
defined, it is sensitive to systematic error and selection that 
may have tended to exclude hotter (or colder) clusters. 

3.2 The X— ray Luminosity Function 

The second X-ray observable used in this paper is the clus- 
ter X-ray Luminosity function. For simplicity we adopt a 
Schechter function parameterization 



dn _ ( L x 

dT x V Li 



l: 



(3.2) 



The local XLF used in this paper is the one given by 
Ebeling et al. Jl997| ) based on the ROSAT Brightest Clus- 
ter Sample (BCS), which contains clusters with z < 0.3, 
selected and flux-limited at X-ray wavelengths. We adopt 
their best fit parameters to equation for the 0.5 — 2.0 keV 
band, namely [A = 33.2 (10 -8 Mpc-^AL" 1 ); a = 1.85; L* = 
5.7 (10 44 ergs" 1 )] and [A = 49.5 (10~ 8 Mpc^A/," 1 ); a = 
1.82; L* = 10.7 (10 44 ergs" 1 )] for the 0.3 - 3.5 keV band. 
The functions are plotted as solid lines in Figure |2| Since 
they find no significant evolution (< 1.8cr) in their sample, 
we assume this to represent the XLF at the present day. 

In order to place constraints on the observed evolution 
of the XLF, we have used the distant cluster redshift dis- 
tributions and luminosity functions that are presently avail- 
able, summarised in Table ^[ Figure ^| shows the positions of 
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□ Southern SHARC, 0.30<z<0.70 
X ROSAT DCS, 0.50<z<0.85 
A WARPS 0.30<z<0.70 



.0 -0.5 0.0 0.5 1.0 

log,, L, (0.5-2.0 keV) [10" erg s~'] 




-0.5 0.0 0.5 1.0 

log,, L, (0.3-3.5 keV) [10 4< erg s - '] 



.<0.3|ROSAT BCS 




-0.5 0.0 0.5 

log,, L„ (0.3-3.5 keV) [10 4 ' 



erg s ] 



Figure 2. The X-ray luminosity functions used in this paper 
compiled from data observed using the ROSAT and Einstein 
satellites. The solid line in all three plots is the best-fit Schechter 
function to the ROSAT BCS z < 0.3 sample as determined by 
Ebeling et al. (1997). Plot (a) illustrates the high-redshift non- 
parametric XLF data from ROSAT surveys, evaluated in the 
0.5 — 2.0 keV pass band. The square points are from the Southern 
SHARC survey (Burke et al. 1997), the crosses from the RDCS 
survey (Rosati et al. 1998, private communication) and the trian- 
gles from the WARPS survey (Jones et al. 1998, private commu- 
nication). Plot (b) illustrates the EMSS cluster X-ray luminosity 
function from the work of Henry et al. (1992) where plot (c) is 
the cumulative XLF point from the EMSS distant cluster sample 
(JV(> 5 X lO^crgs- 1 ) = 1.1 ± 0.43 X 10~ 8 Mpc- 3 ; Luppino & 
Gioia 1995), both evaluated in the 0.3 — 3.5 keV band. The vertical 
dashed line marks the position of L*. Data is presented assuming 
g = 0.5 (fi = 1). 



the high-redshift data, evaluated for qo = 0.5(f2 = 1). Fig- 
ure |j|(a) illustrates the binned, non-parametric XLF's that 
have been compiled from ROSAT data (SHARC , WARPS 
and RDCS ) , with luminosities evaluated in the 0.5—2.0 keV 
band. None of these samples show significant evolution of the 
cluster XLF, out to typically z ~ 0.5. Figure ^(fe) shows the 
EMSS 1 XLF in the 0.3 - 3.5 keV band. Clearly, this hints 
at negative evolution (i.e. a lower space density of clusters 



of given luminosity at higher redshift), although one must 
note that higher luminosity bins usually have greater me- 
dian redshift values and hence applying one redshift to the 
whole sample leads to an overestimation of the evolution. 
Finally, Figure ^(c) shows the EMSS 2 cumulative point, 
also in the 0.3 — 3.5 keV band. Since this point is signifi- 
cantly lower than the BCS XLF, this should provide a tight 
constraint on n and e. 

We use a likelihood method to compare the data-points 
with each model prediction defined by (n, e), for the assumed 
cosmology. Ideally, we would constrain our model parame- 
ters on the basis of individual clusters using a global maxi- 
mum likelihood approach. Starting from the local XLF, we 
could assign individual likelihood probabilities to each of 
the observed clusters (and non-detections) by combining the 
model X-ray luminosity function at the appropriate redshift 
with the selection function defined by each of the surveys. 
This approach would avoid all problems related to redshift 
binning of the available data, and allow cosmological correc- 
tions to be consistently applied. Unfortunately, the detailed 
survey selection functions are not generally available to us, 
so we must adopt an approximate approach. Since little evo- 
lution is observed, binning the data in redshift is not likely 
to result in significant bias. 

The constraints placed on the model parameters (n, e) 
from evolution of the various high-redshift XLF's (relative 
to the BCS XLF) for a given cosmology were generated as 
follows. Firstly, we converted the high-redshift XLF data 
points to agree with the assumed cosmology. In Qq < 1 
cosmologies the predicted space density of clusters decreases 
because the luminosity limit of the survey is brighter (scaling 
as 1/Knax, where Knax is the volume limit of the clusters) , 
hence the luminosities of the clusters themselves are higher 
(scaling as d 2 L , where <1l is the luminosity distance). We scale 
the space densities assuming that all clusters are at the flux 
limit of the survey, however the correction corresponding to 
clusters being a factor of two brighter is ~ 1%. We then 
assumed that the number of clusters in the i* bin (N c ius,i) 
were sampled from a Poisson distribution, with the expected 
number being a function of the true values (n T ,e T ). We 
calculated the C-statistic for the range of model (n, e) such 
that —2.5 < n < 1.0 and —6 < e < 4, where C is defined as 



Cxlf 



N 



In Pi (Nci us ,i;Ei(n,e)), 



(3.3) 



N is the number of bins in the high-redshift XLF and Pi 
is the probability of observing N c i uat i clusters given Ei, the 
expected number calculated from the model, at the median 
redshift of the data point. Confidence regions were then cal- 
culated by differencing C with the minimum value (i.e. the 
most probable set, (no,eo)), AC = C — C m in(no, eo)- If (n,e) 
were independent parameters, AC would be distributed as 



X with two degrees of freedom (Cash 1979), but we do not 



assume here that this is the case. To circumvent this prob- 
lem, we generated a large number of monte-carlo realisations 
of each XLF data-set, for all the cosmologies studied. Each 
realisation involved generating a set of XLF points at the 
same luminosities as the real data-set, drawn from Poisson 
distributions with their mean values set to the numbers of 
clusters predicted by the best-fit model. We then took the 
C distribution for each realisation and calculated the dif- 
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ference in C values between this distribution and the one 
produced by the best-fit model, at the fixed point (no,eo)- 
This method allows us to build up a likelihood distribution 
of AC values that resembles the true distribution. Confi- 
dence levels were chosen such that they enclosed 68% and 
95% of the total number. 



3.3 Results 

For each of the fiducial values of fio (in both open and flat 
Universes), we determined the likelihood distribution as a 
function of n and e. Initially, we analysed the data-sets in- 
dependently, however the results show similar trends, albeit 
less constrained than the overall likelihood distribution, gen- 
erated from combining all of the XLF data. We therefore 
focus our discussion using the combined XLF data-set - the 
result of this is shown in Figure ^. Clearly, this illustrates 
how the constrained area of parameter space changes with 
the background cosmology. This is the result of an interplay 
between two contributing factors, namely the model scaling 
relations used to fit the high redshift data and the positions 
of the data-points themselves, relative to the local XLF and 
L x — T relation. 

Constraints from the evolution of the L x — T normal- 
isation show a weak degeneracy between n and e for high 
values of fio. As fio —> 0, the range of e becomes indepen- 
dent of n (centred on e — 0). The range of n is uncon- 
strained for all values of fio. The reason for the degener- 
acy is straightforward. The data imply no net evolution of 
the L x — T normalisation for an empty Universe (fio = 0). 
Hence, increasing the value of fio pushes the high redshift 
sample to lower luminosities, prompting positive evolution 
in temperature for a fixed value of L x . In an fi = 1 Universe 
with no entropy evolution (e = 0) and a power spectrum 
with a negative slope, the model over-predicts the cluster 
luminosities, hence positive entropy evolution is required to 
match the data. As n increases, the over-prediction becomes 
less severe (hence e decreases). For fio < 1 cosmologies, the 
amount of evolution required becomes progressively weaker 
as fio decreases. Since the growth of structure is also slower, 
the data becomes more and more consistent with e = for 
all values of n. 

The results from constraining the evolution of the com- 
bined XLF data-sets show a degeneracy in the opposite 
sense: higher values of n demand higher values of e. This 
relation can be clearly explained from the fi = 1 plot alone. 
Again, the model scales the XLF in such a way that for 
e = 0, the result is an over-prediction of the cluster lumi- 
nosities. However, this effect becomes less pronounced as n 
decreases, leading to lower values of e. Eventually, the local 
luminosity function is shifted so much that the curvature 
becomes important, leading to an under-prediction of the 
abundance of clusters at the faint-end, where most of the 
data-points lie. This effect places a lower limit on the value 
of n. For lower values of fio, the data-points are pushed to 
brighter luminosities, as well as lower values of dn/dL x . Al- 
though the faint-end data are consistent with no evolution 
for fi = 1, the re-scaling forces the points below the local 
XLF for lower values of fio. The brighter data (in particular 
the EMSS 2 point) is less sensitive to the change in cosmol- 
ogy, due to the steeper slope of the XLF. Coupled with the 



fact that the evolution of structure is weaker, an adequate 
fit to the data is allowed for more negative n. 

In comparison to the open models, the flat (fio+Ao = 1) 
models show slight variations in the constrained region of 
parameter space. The reason for the difference is two-fold. 
Firstly, the data-points themselves are scaled differently be- 
tween cosmologies than their open counterparts. This is be- 
cause the distance to a fixed redshift is larger for flat models 
(hence also increasing the volume element). This has the ef- 
fect akin to having an open model with a lower value of fio. 
Secondly, the growth of structure in a flat model is more 
rapid than the corresponding open case, which affects the 
model scaling relations not unlike an open model with a 
higher effective fio. 

Combining the constraints from the L x — T and XLF 
measurements, we are unable to determine the value of fio 
(with or without Ao) from the data alone: there is always 
a region of parameter space consistent with both observa- 
tions. The same conclusion applies to all of the individual 
XLF data-sets. Specifically, the faint -end data-points from 
the R OSA T surveys are unable to set tight constraints; how- 
ever they are able to set a lower limit on the slope of the 
power spectrum that is sensitive to the value of fio- The 
EMSS samples, particularly the EMSS 2 cumulative point, 
start to provide tighter limits on n and e, since these results 
have started to probe the more sensitive luminosity range, 
brighter than L*. Therefore, future surveys will be more ef- 
ficient at extracting information on the evolution of clusters 
by covering larger areas of sky, probing luminosities around 
and beyond L x . 



3.4 Breaking the n— fio degeneracy 

It is evident that the model is unable to determine a value 
of fio from the XLF and L x — T normalisation constraints 
alone - a feasible range of n and e can always be found for the 
particular scenario. Noting that the constrained range of n 
is sensitive to fio, combining an independent determination 
of n with the X-ray data will be able to place a significant 
constraint on the value of fio. There are two ways in which 
this can be achieved: a theoretical prediction of the primor- 
dial power spectrum that gives the (cosmology-dependent) 
slope on cluster scales or a direct measurement from redshift 
surveys. 

For the theoretical case, we need to obtain the local 
slope of the power spectrum on cluster scales, n e ff such that 



dlnP 



"-off 



dlnfc 



(3.4) 



Although there has been feverous debate recently, the cur- 
rently most popular cosmogony is that dominated by cold 
dark matter (CDM), with n — > 1 for small k, turning over 
to n ~ —3 on small scales. We have adopted the transfer 
function given by Bardeen et al. (1986). The turnover scale 
is determined by the shape parameter, F, which is in turn, 
determined by the cosmological model: 

r = fioh exp 



(3.5) 



I 



ugiyama 1995 ) . We take the value fib — 0.013 h 2 from nu - 
cleosynthesis constraints (Copi, Schramm & Turner 1995) 



No contribution enters from the cosmological constant, due 



© 0000 RAS, MNRAS 000, 000-000 



8 Scott T. Kay and Richard G. Bower 



4 






2 




. - 









UJ 






-2 


7 1 




-4 


1 1 




-6 


: , 1 





-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 

n 



-2.5 -2.0 



^ \ ' 




- 1 








: 1 
: 1 
: 1 


n =o.5 : 



1.5 -1.0 -0.5 0.0 0.5 1.0 
n 



4 


X, 1 




2 


: "r 







'- 1 




UJ 


""7* 




-2 


: "\ ' 




-4 


: 1 


O o =0.5 ; A =0.5 I 


-6 


: 1 , 





-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 



4 

2 




r _ : 


4 

2 




1--.. 
































UJ 










-2 






-2 












i 








\--- 






-4 




i 


n =o.3 : 


-4 




\ 




O o =0.3 ; A =0.7 : 


-6 




i , 




-6 




1 , , . 







-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 

n 



-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 

n 




n„=o.i 



-4; [ 

-fi t. . . I, ■ . I 1 

-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 
n 




-2.5 -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 



Figure 3. Constraints on the slope of the power spectrum (n) and the entropy evolution parameter (e) for seven different cosmological 
scenarios. Solid contours are the 68% and 95% confidence regions from the evolution of the X— ray luminosity function, using the combined 
data from all of the surveys. Dotted contours represent the 68% and 95% confidence regions for evolution of the cluster luminosity- 
temperature relation. Shaded regions are for parameters consistent with the evolution within 68% li mits (of both XLF and L x — T). The 
vertical solid line is the determination of n from the CDM based theoretical calculation in Section 3.4 and the vertical dashed line is a 
measured value of n from the APM rich cluster sample. 
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to the fact that it contributed negligibly to the energy den- 
sity at the last scattering epoch. To estimate fc e ff , we assume 
that the virial mass of a typical rich cluster in the local Uni- 
verse (e.g. the Coma cluster) is M ~ 5 x lO 14 h _1 M . If 
the cluster formed by the top-hat collapse of a spherically 
overdense region, the comoving radius required to contain 
this mass (assu ming present background density) would be 
R ~ 8 h" 1 Mpc (|Eke, Cole & Frenk 19961). This is only true 



for a critical density Universe but can easily be modified for 

_ i_ 

lower densities, since R oc fi 3 . The value of R can then be 
converted to an effective scale o n the linear power spect rum 
by using the following formula ( Peacock fc Dodds 1994 ) : 



;|(n+l)]! 



R ' 



(3.6) 



which should iteratively con verge on the values (n e ff,fc e ff) 



with the use of equation 3.4. The calculated values of n are 
represented as solid lines in Figure ^. Note that the value 
of n becomes more negative with decreasing values of Qo- 
Although the effective scale of galaxy clusters is larger in 
a low mass-density Universe, the change in shape of the 
CDM power spectrum itself (parametrised by F) dominates 
the shift in the value of n. As Qo decreases, the position 
of the turnover moves to larger scales, leading to a more 
negative slope on cluster scales. Figure |^ shows the result 
of this calculation, plotted as a solid vertical line. Analysing 
the likelihood distributions for a continuous range of Qo, we 
find that the CDM prediction is only consistent with the 
X-ray data for Q < 0.55 (A = 0) and 0.1 < Qo < 0.7 
(Qo + Ao = 1), to 95% confidence. We can place a lower 
limit on Qo in the flat model because the predicted value 
moves to more negative n faster than the XLF contours, as 
Qo decreases. 

Another approach is to constrain n directly by measur- 
ing the power spectrum on rich cluster scales from wide 
area sk 
ton 19! 



Dal- 



surveys. One example (Tadros, Efstathiou 
comes from the APM cluster sample, which has 
a median redshift (z) — 0.09, measuring a slope, tiapm = 
— 1.6 ± 0.3. We currently assume that the underlying power 
spectrum has negligible curvature on the range of scales we 
are studying. Ideally, we would separate the population into 
redshift bins and scale each respective population relative 
to the preceeding one, using the slope of the spectrum at 
that epoch. Imminent surveys such as the Sloan Digital Sky 
Survey and the 2dF Galaxy Redshift Survey should provide 
the necessary data to constrain the slope at higher redshifts. 
Furthermore, the APM result assumes that the bias is linear 
and therefore scale-invariant on scales of the mean inter- 
cluster separation. 

To constrain the range of Qo, we took the 95% limits 
placed on n and e from the evolution of the XLF and L x — T 
normalisation and subsequently fixed the value of n to —1.6. 
For the XLF data, we found that the model fits the data 
more and more adequately as Qo decreases (for the range of 
n considered): this is shown in the plot by the shift in the 
contours. Hence, this prevents us from placing a lower limit 
on Qo- In the flat models the geometry weakly alters the 
shape of the likelihood distribution but the overall trend 
is the same. Since both constraints produce contours that 



Survey 


Ho, 


u(A = 0) 


f2 , u (A 


o = 1 - Qo) 


SHARC 


1.0 


[0.6, > 1.0] 


1.0 


[0.5, > 1.0] 


WARPS 


1.0 


[0.75, > 1.0] 


1.0 


[0.7, > 1.0] 


RDCS 


1.0 


[0.65, > 1.0] 


1.0 


[0.6, > 1.0] 


EMSS 1 


1.0 


[0.7, > 1.0] 


1.0 


[0.6, > 1.0] 


EMSS 2 


0.7 


[0.4, > 1.0] 


0.65 


[0.3, > 1.0] 


COMBINED DATA 


0.4 


[0.25,0.6] 


0.4 


[0.2,0.65] 



Table 3. Upper limits on the value of f2o (with and without 
Ao), using the constraint that n = —1.6, from the APM survey 
(Tadros, Efstathiou & Dalton, 1997). Results from both individ- 
ual and combined XLF data— sets are presented. Limits are set 
using the 95% confidence limits consistent with both the XLF 
and L x — T constraints. Ranges in square brackets are for quoted 
uncertainties in n (±0.3) 



become more aligned in lower Qo cosmologies, we are unable 
to make further conclusions regarding the lower limit. 

What the data does allow us to do is place an upper 
limit on the value of Qo- The results are given in Table ^, 
displaying the values that are consistent with the 95% con- 
fidence limits of the XLF and L x — T relations. Results for 
both the individual and combined XLF data-sets are pre- 
sented. The ranges in brackets were determined by varying 
the value of n by ±0.3. It is evident that the value of Qo is 
unconstrained using the APM value, for all of the individual 
data-sets except the EMSS 2 survey. The flat models give 
similar results to their open counterparts, with the range in 
the upper limits being slightly larger. Combining the XLF 
surveys sets a tighter constraint: Qo < 0.6 (or equivalently 
Qo < 0.65 for Qo + Ao = 1). 

3.5 Evolution of Cluster Core Sizes 

Finally, we consider another possible constraint that could 
be readily measured from X-ray data: the evolution of clus- 



ter core-radii. Using the model we developed in section 2.3 
we predict the characteristic core-radius of a cluster popu- 
lation to scale as: 



(3.7) 



b = 



d 



a - i V. ' /3 

20' 

2^ _1, 



3 V 2/3 



Even though the core radius depends on both n and e, 
it's evolution is determined much more sensitively by the 
changes in the central entropy of the gas rather than the 
amount of structural evolution. Figure illustrates the re- 
gions of parameter space consistent with various amounts of 
core radius evolution, out to a redshift z — 0.5. The back- 
ground cosmology and spectral index have only a weak influ- 
ence on the evolution of r c , hence observations should quite 
easily place a constraint on the value of e - the two extremes 
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Figure 4. The region of parameter space covered by evolution 
of the cluster core radius for various models. The shaded regions 
depict the range of cosmologies, from Q = 1 (with the great- 
est amount of structure evolution) to Q = 0.1, Ao = 0.0 (with 
the least amount of structure evolution). Each particular area il- 
lustrates the different amount of evolution in the redshift range 
z = [0.0,0.5], with the labels indicating the value Alogr c . 



of the shaded regions represent scaling in the core by a fac- 
tor of ~ 0.3 and 3 respectively. Core sizes that shrink with 
redshift demand negative values of e (i.e. the core entropy of 
clusters decreases with redshift) whereas core sizes that grow 
with redshift demand values of e that are weakly negative 
or higher. This approach is useful in the sense that we can 
directly measure the rate of entropy evolution in the core 
and test it's consistency with the other methods presented 
above. 

Analysis involving looking at changes in the cluster core 
size has now started to appear in the literature. Results are 
currently suggesting no si gnificant evolution in t he core size, 
at least out to z ~ 0.5 (Vikhlinin et al. 1998). Given the 
uncertainties, our results are consistent with this, however 
positive evolution in the cluster core radii (e > 0) would be 
preferred. At present, the sample sizes are small: a more 
comprehensive measurement of core radius evolution will 
serve as an important test of the validity of this model. 



4 DISCUSSION 

At this stage, we go on to focus on several issues that are 
worthy points for discussion. Firstly, our results are based 
on several important assumptions made when modelling the 
X-ray evolution of clusters, so we discuss our choices and 
their robustness. Secondly, particularly due to the recent 
publication of ROSAT and ASCA data, several results con- 
straining Qo using the evolution of X-ray clusters have now 
appeared in the literature. We therefore compare our results 



with those found by other groups. Finally, we address what 
we think is an important strategy for observations (partic- 
ularly with the imminent launch of the XMM satellite, and 
pin-point the areas that we think will maximise future suc- 
cess in this field. 



4.1 Model Assumptions 

Although general and based on physically motivated scal- 
ing relations, there are several key approximations in the 
model used to arrive at this conclusion. One example is the 
(3 profile: the results presented in §3 assume that the X-ray 
emission profiles are well fit when using the j3 — model, where 
the density of the intracluster gas tends to a constant value 
at small r (the core) and the isothermal result at large r, for 
a value of 2/3. Individual measurements of (3 do not neces- 
sarily give this value, although we assume that the average 
value over the population at a given epoch is well repre- 
sented by ((3) = 2/3. However, it is important to investigate 
the tolerance of the model when varying the value of /3. Fig- 
ure ^ illustrates the effect on the constrained range of n and 
e that variations in (3 induces. Using the SHARC data as 
an example, we have taken two values of (3 : 0.6 and 0.7, 
defining what we think are the expected limits of any fluc- 
tuation. Two cosmological models have been chosen : Q — 1 
and Qo = 0.1, which represent the most and least amount 
of structural evolution respectively. Changing f3 can signifi- 
cantly change the range of e values but has relatively little 
effect on the values of n. This, coupled with the work on 
the core radius discussed in the last section, reinforces the 
idea that both n and e contribute separately to the evolu- 
tion of clusters: the entropy evolution is sensitive to the way 
in which the gas is distributed, with /3 and r c defining the 
shape of the gas density profile, while the mass evolution 
is tied in with the background density. Hence, the depen- 
dence of Qo on (3 is weak, although certain models may be 
ruled out on the basis that the range of e may be implau- 
sible. This is evident in the figure, where the combination 
((3 — 0.6, Qo — OA) pushes allowed values of e into the limit 
set by the maximum cooling calculation (e = 2.0) given in 
Paper I. 

With the Launch of the AXAF satellite, it will become 
possible to study the emission profiles of distant clusters in 
more detail. These studies may show that the cluster profile 
evolves with redshift. For example, simulations of c luste r 
evolution with pre-heated gas by Mohr & Evrard (1997) 
suggest that the entropy injected reduces the overall slope 
of the gas density profile rather than just that of the profile 
of the very centre. This effect can be parameterised by a 
much slower role over near the core radius but the scaling of 
the solution remains the same, since the X-ray emission still 
scales as p^,r^T a . Thus, as Figure |B| shows, the j3 profile and 
the e parameter are interconnected - while observations of 
the emission profiles of high redshift clusters would cause us 
to revise our physical interpretation of a particular e value, 
and re-examine the balance of heating and cooling within 
the cluster, it would not invalidate the phenomenological 
description provided by the model. 

We may apply the same argument to take into ac- 
count the non-is othermal temperature of structure of clus- 



ters suggested by 



Markevitch et al. 1998 



The model only re- 



quires that the temperature structure is normalised by the 
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Figure 5. An example of the effect of varying /3 using the SHARC XLF data. The left column is for Q = 1 (the most mass evolution) and 
the right column Qo = 0-1 (the least mass evolution). The top and bottom rows for the expected limits on f3 (0.6 and 0.7 respectively). 
The effect on n is insignificant, although raising /3 lowers the range of e. 



virial temperature of the dark matter halo. The tempera- 
ture structure may affect the rate at which the cluster core 
evolves for a given entropy change, but this should again 
taken into account when the e parameter is interpreted. 

The model has been developed under the assumption 
that clusters of galaxies are in dynamical equilibrium. Clus- 
ters that are seen during major mergers may have X-ray 
luminosities and average temperatures that are far from 
their hydrostatic values. Nevertheless, the description we 
have developed applies to the cluster population as a whole, 
and the properties of individual clusters may differ wildly 
from the average. In this sense, the effects of mergers are al- 
ready incorporated through the entropy evolution described 
by the model, and the disruption of hydrostatic equilibrium 
is mimicked by an increase in core entropy. The cluster tem- 
peratures are also susceptible to these departures, however, 
since they are weighted by the luminosity. Ideally, we should 



use clu ster temperatures measured from the outer regions 



of clusll oro whore the gao temperature bettor roflocto the 
virial temperature of the gravitational potential. These are 
not available for distant clusters at present, but our results 
will still be valid so long as the evolution of the luminosity 
weighted temperature reflects the evolution of the cluster 
potential. This is a problem common to all methods based 
on the X-ray temperatures of clusters. 

The slope of the luminosity-temperature relation does 
not explicitly enter our determination for the evolution of 
clusters. All that is required in our model is the evolution 
of the normalisation. We have implicitly assumed that the 
slope remains fixed at its present-day value. This assertion, 



which corresponds to the assumption that the entropy evo- 
lution parameter (e) is independent of cluster temperature, 
is unavoidable with the limited temperature data that is 
available for high redshift clusters. Consequently, it is impor- 
tant that the evolution of the normalisation is determined 
by using consistent slopes for high and low redshift data, 
and preferably using high redshift clusters representative of 
the clusters used to determine the luminosity function evo- 
lution. At present, temperature measurements are available 
only for the brightest clusters at high redshift. This situation 
will hopefully improve with the launch of the XMM satel- 
lite (Lumb et al. 1996); its large collecting area is ideally 
matched to the determination of temperatures for the lower 
luminosity, distant clusters. The slope of nearby (z < 0.05) 
cluste rs has been most recently determined by Markevitch 
( 1998| ), using 7? OS^T luminosities and ASCA temperatures, 
measuring a value, A ~ 2.65. However, we make use of 



the olde r sample based on Einstein MPC data (David et 



al. 1992 ) for two reasons. Firstly, constraints on the evo- 
lution of t he L x — T relation are provided by Mushotzky 



& Scharf (1997), which uses this (local) sample to contrast 
with their higher redshift ((2) ~ 0.3) ASCA data. Also, the 
slope supplied by Markevitch has been calculated for clus- 
ters with cooling flows removed, in order to directly compare 
with models that exclude non-radiative components in the 
plasma evolution. Since our model is based on evolution of 
the central entropy, we cannot exclude the cooling flow clus- 
ters, which probably contribute to most, if not all of the 
negative entropy evolution in cluster populations. We argue 
that the dominant cause for the intrinsic scatter in the re- 
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lation is due to each cluster having it's own particular core 
entropy, coming from their individual formation histories. 

The considerations we discussed above related to the 
distribution of gas within the dark matter potential. In a 
low density universe, it is necessary to explicitly account 
for the distinction between the background density of the 
Universe when the cluster is formed, and it's value when 
the cluster is observed. We have related these two epochs by 
assuming that the mass of the cluster grows by a factor of 
two (Kitayama & Suto, 1996). However, for fio > 0.1, this 
factor produces only a small change in the halo evolution. 

4.2 Comparison with Other Methods 

A number of other papers have recently appeared in the lit- 
erature dealing with the constraints on fio from the X-ray 
evolution of clusters. One of the most popular methods uses 
the evolution of the cluster abundanc es (e.g. [Henry 1997j 
Eke et al. 1998|, IViana & Liddle 199S|, iMathiesen & Evrard 



ical justification. Reichart e t al. a dopt the same parameter- 
ization as Evrard & Henry (1991), notably L x oc T p (l + z) s ; 



199Sj , pukbir fc Blanchard 1997| , |R.cichart ct al. 1998| ). To 
extract the value of fio, an estimator for the cluster virial 
mass is required. This can be achieved by calculating the 
cluster temperature function at any given epoch and con- 
verting this to mass using the Press-Schechter formalism 
with the assumption that the cluster haloes form via spheri- 
cal infall. The result from this method is still slightly unclear 
- at present the majority of authors find values of fio typi- 
cally in the region 0.3—0.5 with or without A, alt hough some 



auth o rs find higher values includin g fi = 1 (e.g. Sadat et al 
199S , Blanchard & Bartlett 1998). An extensive discussion 



on the exist ence of this clash in results can be found in Eke 
et al. (1998). In particular, they highlight the discrepancy 



in the abundance of low-redshift clusters, particularly the 
high temperature end, where the measured value is uncer- 
tain. Those favouring high fio assume a higher abundance 
of clusters in this region than those favouring low values of 
fio. 

Our work differs from these treatments in that it does 
not assume a detailed model for the abundance of clusters 
or for their temperature or X-ray luminosity distribution. 
Instead, we have emphasised the scaling relations that must 
relate the properties of clusters at one epoch with those at 
another. By introducing an additional parameter that de- 
scribes the evolution of the core gas entropy, our approach 
allows us to inter-compare clusters at different epochs with 
minimal further assumptions. Comparing our results with 
those obtained by other authors illuminates the constraints 
which are general and th ose wh ich are model-specific. As an 
example, Reichart et al. (1998) also look at luminosity func- 



tion evolution by adopting the Press-Schechter approach for 
the evolution of the mass function , but use the L x — T re- 
lation to calibrate their mass to X-ray luminosity relation. 
In addition, their model differs from ours in that it involves 
the slope of the L x — T relation explicitly, resulting in a 
purely empirical calibration of the luminosity-mass conver- 
sion. By contrast, our approach uses the weak self-similarity 
principle (plus the entropy evolution) to scale the properties 
of clusters between epochs. As discussed extensively in Pa- 
per I, this is completely separate from the scaling relations 
between clusters of different mass at the same epoch. Try- 
ing to equate the two requires that we adopt a Strong Self- 
Similarity principle; something for which there is little phys- 



assuming (3 = 2/3, we find for our model that p — 13/6 
and s — (13 — 3e)/4. However, as we have emphasised, the 
way in which we constrain our parameters (particularly e in 
this context) is fundamentally different, and we should not 
expect to obtain identical results. 

We choose to co ncentrate our comparison with the work 
of Eke et al. (1998). These authors derive a value of fio 



from the evolution of the cluster temperature function, using 
emission-wei ghted temperatures obtained from ASCA data 
(Henry 1997). Since the cluster sample is X-ray flux limited, 
they have to correct for incompleteness in their estimation 
of the temperature function. This is done by assuming a 
non-evolving L x — T relation (i.e. fixed slope, normalisation 
and scatter) . The advantage of this method is that it infers 
the mass distribution in a more direct manner than our ap- 
proach, hence their work places a strong constraint on fio: 
fio = 0.43 ± 0.25 for A = and fi = 0.36 ± 0.25 with 
a non-zero cosmological constant. This was achieved with 
a substantially smaller data-set. Furthermore, the authors 
were able to internally constrain the shape of the power spec- 
trum. This is possible because they directly link the form of 
the temperature function to the mass function through the 
Press-Schechter model (assuming a CDM power spectrum). 
They arrive at a value of V of 0.08 ± 0.07 (r = 0.09 ± 0.08 
for Qo+Ao = 1). 

To compare our results to those of Eke et al., we first 
convert their measured range of P to an equivale nt r ange 
of n, using the same method discussed in section 3.4. We 



then find the range of fio that is both consistent with the 
X-ray data and the calculated range of n. We find that 
fio < 0.45(fi < 0.5 for fi + A = 1) to 95% confidence, 
in agreement with their upper limit. However, we fail to 
place a constraint on the lower limit of fio. This is initially 
surprising, but seems to be due to systematic differences in 
the data. We find that, for a low fio cosmology, the available 
data are consistent with weak negative evolution of the XLF 
and no evolution in the L — T normalisation. By contrast, 
Eke et al.'s temperature function data (corrected for cos- 
mology) shows a decline in temperature function amplitude 
suggesting that a very low value of fio is unacceptable due to 
the deficit of clusters at this epoch. Such dependence on the 
particular data set is worrying, but emphasises the impor- 
tance of using independent data-sets to establish consistent 
conclusions and reduce systematic biases. It is satisfying to 
find that our method finds similar results, at least in deter- 
mining the upper limit. 

4.3 Future Directions 

Since the results presented in this paper were based on scal- 
ing the ROSAT BCS XLF to fit the higher redshift surveys, 
we have assumed that the two populations are well sepa- 
rated in redshift. This is not immediately obvious, since the 
ROSAT BCS sample has clusters out to z ~ 0.3, where the 
other samples start around this value. The saving grace is 
that the higher-redshift surveys have detected lower lumi- 
nosity clusters (L x < L x ), due to the small solid angles of 
the surveys, whereas the BCS clusters in this regime are 
at much lower redshifts. However, this range of luminosities 
places a much weaker constraint on the model than if high 



© 0000 RAS, MNRAS 000, 000-000 



Measuring fi from the Entropy Evolution of Clusters 13 



redshift data were available at the brighter end (L x > LJ). 
Improved coverage will result in a better internal constraint, 
since the effects of evolving cluster number density and clus- 
ter luminosity will be more easily separated. As a pointer for 
the future, what is required is a more statistically sound de- 
scription of the high-redshift XLF from L* and above. As an 
hypothetical example, we assume that the combined data- 
set used in this paper is a good description of the XLF for 
luminosities around and below L*. A measurement is then 
made at L x (0.5 - 2.0 keV) ~ 10 45 ergs" 1 , detecting 10 clus- 
ters at (z) = 0.5 in a solid angle that resulted in a cumula- 
tive abundance in agreement with the BCS cumulative XLF 
(N(> L x ) ~ 0.25 x 10" 8 Mpc" 3 ). This strongly suggests a 
detection of no evolution, particularly for the Q = 1 geome- 
try. Using the APM constraint on n, we find that Qo < 0.2 
at the 99% level. In practice, the XMM satellite is particu- 
larly suited to this task, with its large collecting area and 
high s patial resolution, that will allow one to easily make 



the distinction between clusters and AGN. What is required, 



however, is a wide area survey covering at least 500 square 
degrees, rather than a survey that is particularly deep. In 
the near term, progress will be made through the ROSAT 
North Ecliptic Pole survey - a wide angle surv ey composed 



of s uper posed ROSAT All-Sky Survey strips (Henry et al 
199SpT| 



5 CONCLUSIONS 

In this paper, we have taken the entropy-driven model of 
cluster evolution, detailed in Bower (1997) for an Q = 1 
Universe, and modified it to account for evolution in differ- 
ent cosmological scenarios, specifically open models and flat, 
non-zero A models. This allows us to separate contributions 
made by the hierarchical growth of structure (controlled by 
the slope of the power spectrum, n) and changes in the core 
entropy of the intracluster gas (controlled by the entropy 
evolution parameter, e). 

We then placed constraints on n and e for seven reason- 
able cosmological models, using the current wave of X-ray 
data for the two observables: the X-ray Luminosity Function 
and the Luminosity-Temperature relation. For the L x — T 
relation, we have taken the compiled low and high redshift 
samples from David et al. (1993) and Mushotzky & Scharf 
(1997) respectively and determined the best-fit slope and 
scatter. We subsequently placed 68% and 95% confidence 
limits on the evolution of the temperature normalisation, 
which is slope-dependent. For the XLF, we have used the 
r ecently available ROSAT sa mples, specifically the SHARC 
flBurkc et al. 1997]), WARPS (pones et al. 1998|) and RDCS 



(Rosati et al. 1995) high-redshift, non-parame tric determi- 



nations, as well as the older EMSS samples (Henry et al 



1992, Luppino & Gioia 1995). The amount of evolution in 



these measurements was q uantified by compar ing with the 
local ROSAT BCS XLF (lEbeling et al. 19971). Using the 



X-ray data alone, we find acceptable regions of parameter 
space for every cosmological model considered, to a high 
degree of confidence. Specifically, as the density parameter 
decreases, more negative values of n are allowed to com- 
pensate for the weakening growth rate of clusters. Values of 
constrained e are appropriate to modify the luminosity evo- 
lution to compensate for the measured change in abundance, 



moving weakly from positive values to zero (no evolution) 
for low-density Universes. 

To break the degeneracy between n and Qo, we calcu- 
lated the slope of the power spectrum based on the CDM 
paradigm, by taking the typical mass of a rich cluster and 
calculating the corresponding fluctuation scales that give 
rise to such objects. We find that the CDM calculation sets 
a limit of fi < 0.55 (0.1 < Q < 0.7 for Q + A = 1). Using 
the value determined by the APM survey (n = —1.6 ± 0.3), 
we conclude that Qo < 0.6 (< 0.65 for Q +Ao = 1). All lim- 
its were calculated to 95% confidence. In general, we cannot 
determine a lower limit to the value of Qo , since our model 
fits the data much more comfortably at fixed n, for lower 
density Universes. The evolution of the cluster core-radius 
was also considered, where we showed that it depends much 
more sensitively on the entropy evolution of clusters than 
the structural evolution (determined by n and Qp ). Present 
data sugges t no strong evolution out to z ~ 0.5 (Vikhlinin 



et al. 1998), roughly consistent with our results given the 



quoted uncertainties. 

Finally, we discussed the robustness of assumptions we 
have made in the model, particularly the effect of varying 
the mass distribution of the gas (controlled by the value of 
P). Again, this has a much stronger effect on the constrained 
range of e than the rate of structural evolution, and hence 
only weakly affects the determination of Qq. We then com- 
pared our result s to that of other groups, particularly with 
Eke et al. (1998). We find that our upper limit is consistent 
(f2o < 0.45, or Q < 0.5 for Q + A = 1 at 95% confidence) 
although we cannot place a lower limit on Qo- 

In order to provide tighter limits on the value of Qo, it 
will be essential to obtain an accurate measurement of the 
bright end of the XLF out to redshifts at least comparable 
to present distant cluster samples, ideally from a wide an- 
gle survey. This is something we hope that the forthcoming 
XMM mission can provide. 
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